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Abstract 

Mounting observational data confirm that about 73% of the energy density con- 
sists of dark energy which is responsible for the current accelerated expansion of 
the Universe. We present observational evidences and dark energy projects. We 
then review various theoretical ideas that have been proposed to explain the origin 
of dark energy; they contain the cosmological constant, modified matter models, 
modified gravity models and the inhomogeneous model. The cosmological constant 
suffers from two major problems: one regarding fine-tuning and the other regarding 
coincidence. To solve them there arose modified matter models such as quintessence, 
k-essence, coupled dark energy, and unified dark energy. We compare those mod- 
els by presenting attractive aspects, new rising problems and possible solutions. 
Furthermore we review modified gravity models that lead to late-time accelerated 
expansion without invoking a new form of dark energy; they contain f(R) gravity 
and the Dvali-Gabadadze-Porrati model. We also discuss observational constraints 
on those models and on future modified gravity theories. Finally we review the 
inhomogeneous Lemaitre-Tolman-Bondi model that drops an assumption of the spa- 
tial homogeneity of the Universe. We also present basics of cosmology and scalar 
field theory, which are useful especially for students and novices to understand dark 
energy models. 
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1 Introduction 



All physical bodies attract each other. Then why do our stars, galaxies and the Universe 
not collapse? Einstein believed the Universe is static, but his Einstein equations implied 
a dynamic Universe. In 1917 he introduced the cosmological constant A, the antigravity 
vacuum energy. He had adjusted the scale of A to gravity so that the Universe stays 
static. In the same year 1917 Willem de Sitter has calculated Einstein's equation for a 
vacuum Universe without matter. This solution drives inflation of the Universe. Ironically, 
in Einstein's Universe A was used for a stationary Universe, but in de Sitter's Universe 
A is the reason for inflation. In 1928 Edwin Hubble demonstrated that the Universe is 
expanding and Einstein called the cosmological constant the biggest blunder of his lifetime. 

The discovery of Riess et al. in 1998 [lj and Perlmutter et al. in 1999 [2] changed 
the landscape of a vision of the Universe. They found that distant supernovae at z ~ 0.5 
are Am ~ 0.25 mag, about 25%, fainter than ones expected for a decelerating Universe 
without the cosmological constant. They concluded from the observation of supernovae 
that the Universe is currently accelerating in its expansion. This discovery won them the 
Nobel Prize in physics 2011. 

We call the source of this acceleration dark energy. Not only Type la supernova (SN 
la) observations but also Cosmic Microwave Background (CMB) and Baryon Acoustic 
Oscillations (BAO) substantiate that dark energy captures about 73% of the energy density 
of the current Universe. Dark energy is a major outstanding issue in physics and cosmology 
today. There are a number of useful reviews of dark energy that mainly focused on theory 
[31 IH E], on probes of dark energy [U] and on the cosmological constant [S]. 

The aim of this review is to introduce major theoretical models of dark energy focusing 
on the basic physics of each model for students and non-experts on the subject. Since we 
do not cover all of models in the field, the list of references are limited. The interested 
readers should consult reviews with more complete reference lists, e.g. In the first 
section we present observational evidences (section [T]) . In the following sections we intro- 
duce theoretical approaches to explain the origin of dark energy. In section [2j we review 
the simplest candidate: the cosmological constant A; the cosmological constant suffers 
from two major problems, one regarding fine-tuning and the other one coincidence. To 
solve them there arose modified matter models such as quintessence (section [3]), k-essence 
(section [5]) , coupled dark energy (section [5]), and unified models of dark energy and dark 
matter (section [6]). We compare those models by presenting attractive aspects, new ris- 
ing problems and possible solutions. Furthermore, we review modified gravity models 
which lead to late-time accelerated expansion without invoking a new form of dark energy. 
They include f(R) gravity (section [7]) and the DGP model (section [8]) . We also discuss 
observational constraints on these models. In section [9] we review the inhomogeneous 
Lemaitre-Tolman-Bondi model that drops an assumption of the spatial homogeneity of 
the Universe. Lastly we summarize and give a brief comparison of those candidates. In 
the appendix we present basics of cosmology including the homogeneous FLRW model 
and scalar field theory, which are useful to understand dark energy models. The sections 
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are connected with each other but independent enough for readers to look directly into 
their interested sections without referring to the preceding sections. 

We use the metric signature (+, — , — , — ) and follow mostly the conventions of Mukhanov's 
book [9]. We employ natural units c = h = kg = 1 while retain the gravitational constant 
G = l/M 2 , where M pi = 2.177 x 10~ 5 g = 1.416 x 10 32 K = 1.221 x 10 19 GeV is the Planck 
mass. 

1.1 Observational evidences on dark energy 

The recently dominating dark energy is supported by many independent observations, such 
as SN la [ID] . CMB [11] and BAO JT2] . Here we present how these could be evidences for 
the presence of dark energy. 

Type la supernova (SN la) 

As we mentioned above, SN la observations are one significant evidence for the cosmic 
acceleration. When a white dwarf reaches 1.4M , gravity becomes bigger than the Fermi 
degeneracy pressure and drives the star into explosion. Since SN la is made in the same 
process, we assume that SN la always have the same luminosity and it can be used as 
a standard candle. Measuring the apparent magnitude, we can estimate the luminosity 
distance. A supernova is as bright as a whole galaxy, thus surveying SN la is very powerful 
and a direct way to study the distant Universe. 

But since the detection was first published some scientists had doubt that the distant 
supernovae could appear fainter just due to extinction. Absorption by intervening dust 
could also lead to characteristic reddening. Nowadays thanks to Hubble Space Telescope 
(HST) [T3] we have high quality light curves up to z ~ 1.8 and big amount of other distant 
SN data from ground based observations. We are thus more convinced that this was not 
due to extinction by dust but due to real cosmic acceleration. 

Cosmic Microwave Background (CMB) 

Not only SN la but also CMB observations strongly indicate the presence of dark energy. 
The CMB is a snap shot of the Universe before the cosmic structure developed and has 
temperature anisotropics which are influenced by dark energy. This snap shot tells us the 
cosmic history from the photon decoupling epoch to the present. 

The angular power spectrum of CMB temperature anisotropics measured by Wilkinson 
Microwave Anisotropy Probe (WMAP) [HJ QSJ [161 [17] implies that it is dominated by 
acoustic peaks arising from gravity-driven sound waves in the photon-baryon fluid. And 
then the positions of these acoustic peaks are shifted by cosmic expansion. Thus the 
positions and amplitudes of acoustic peaks contain important cosmic information. 

Baryons were strongly coupled to photons before the decoupling epoch. After recom- 
bination (z r ~ 1100) baryons got free from Compton drag of photons and stayed at a 
fixed radius, i.e., sound horizon. It determines the first acoustic peak in CMB anisotropy 
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that we can observe directly. Since the size of the sound horizon is predicted from theory, 
l s ~ iJ _1 (z r ), this sound horizon serves as a standard ruler by measuring the angular scale 
of the first acoustic peak. It is analogous to how SN la can be used as a standard candle 
but is completely independent from the SN la technique. 

Galaxy clusters also leave an imprint on the CMB, the so-called Sunyaev-Zel'dovich 
(SZ) effect |18j . The SZ effect stands for the distortion of the CMB spectrum through 
inverse Compton scatterings due to the collisions of the CMB photons and the high energy 
electrons in galaxy clusters. Measuring this distortion of the CMB, we can estimate 
masses of the clusters. Since the SZ effect is a scattering effect, its magnitude is redshift- 
independent; very distant clusters are as easy to detect as nearby clusters. When combined 
with accurate redshifts and mass estimates for the clusters, e.g. X-ray observation, the SZ 
effect plays a role of a standard ruler [19] and can be used as a distance indicator. This 
is the same method as comparing the absolute magnitude with the apparent one. The SZ 
effect may also provide interesting constraints on the dark energy. 

The combination of CMB and SN la observations indicates that our Universe is spa- 
tially flat and the matter component is about one quarter of the critical density, fi m ~ 0.3. 
This implies the need for dark energy as a missing component = 1 — — 0.7, which 
is indeed consistent with the value of — 0.7 suggested by SN la surveys [H [2]. 

Dark energy also influences on large angle anisotropies of CMB, the so-called integrated- 
Sachs- Wolfe (ISW) effect [20]. Variation of the gravitational potential during the epoch of 
cosmic acceleration leads to differential gravitational redshifts of photons and large scale 
perturbations. Analyses of spatial correlation of the large scale structure (LSS) of galaxies 
indicate that the Universe is not Einstein-de Sitter Universe (Q m = 1). 

Baryon Acoustic Oscillations (BAO) 

Studying BAO is also a significant investigation for dark energy. As mentioned above, after 
decoupling baryons stay at the distance of the sound horizon whereas dark matter stays 
at the center of overdensity. They attract matter and eventually form galaxies. Therefore 
it is expected that a number of galaxies are separated by the sound horizon. We call it 
BAO signal. It means that by observing the large scale structure (LSS) of galaxies, we 
can measure the sound horizon scales and compare them with theoretical predictions. 

To analyze the BAO signal, we do not need to resolve galaxy magnitudes or galaxy 
images. The only thing to be determined is their 3-dimensional position. Thus in that 
sense, BAO observations are safer from astronomical uncertainties than other dark energy 
probes. 

The Sloan Digital Sky Survey (SDSS) catalog provides us a picture of the distribution 
of galaxies up to z = 0.47 [T2] and search for a BAO signal. To analyze this picture, we use 
a two-point correlation function of galaxies; it gives the probability that one galaxy will 
be found within a given distance of another. The two-point correlation function indeed 
possesses a bump, BAO signal, at a scale lOO/iT 1 Mpc for 0.16 < z < 0.47. Both CMB 
and BAO signals indicate that the sound horizon today is about 150 Mpc. 
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Figure 1: There is strong evidence for the existence of dark energy. Plotted are Q m - Q\ 
(left panel) and Q m - w (right panel) confidence regions constrained from the observations 
of SN la, CMB and BAO. From Ref. [2]]. 



Figure LI shows us the best fit of cosmological parameters according to three inde- 
pendent observations. 



Weak gravitational lensing 

The gravitational bending of light by structure distorts the images of galaxies. Studying 
this gravitational lensing effect, we can map dark matter and its clustering which shows 
us the influence of dark energy in the growth of large-scale structure. 

In rare cases we observe strong lensing which is deflection of light by massive structures 
like clusters and causes multiple images of the same background galaxy. More often we 
observe weak lensing which distorts the shapes, sizes and brightness of galaxies. Since we 
know the shape better than the original size and brightness of galaxy, many studies focus 
mainly on the change on shapes which is called cosmic shear. Weak lensing normally reads 
only 1% change of galaxy shape. Therefore, a large number of galaxies are required to 
detect cosmic shear signal. The shear correlation function can be directly compared with 
the measurements [22J. 

There are dark energy probes other than ones raised above, including cosmic age 
tests from globular cluster and the white dwarfs [23] , galaxy clusters [2U 1251 125] , Hubble 
parameter measurements [23 [2S] and gamma ray bursts [221 EDI EI] • 



5 



1.2 Dark energy projects 

There are many observational projects to detect dark energy using the four dark en- 
ergy probes: SN la, CMB, BAO and Weak Lensing (WL). Among them the Dark En- 
ergy Task Force (DETF) report (32] classifies the dark energy projects into four stages: 
Stage I-completed projects that already released data, Stage II-on-going projects, Stage 
Ill-intermediate-scale, near-future projects, and Stage IV-large-scale, longer-term future 
projects. Table [T] shows us the dark energy projects classified by probes and stages. 



Probes 


SN la 


CMB 


BAO 


WL 


Stage I 


Higher-Z Team [TUl [33], 


COBE [52], 


2dFGRS [56], 


CFHTLS [SO] 




SNLS [33 [351 [36], 


TOCO [53], 


SDSS [57], 


EI] 




ESSENCE [33 [38], 


BOOMERang 


6dFGRS [5H], 






nsf [sg, csp nam], 


[51]. Maxima 


WiggleZ [59] 






LOSS 021 S3J EH], 


[53], WMAP 








SDSS pEi EE], SCP 

[2H mj EH], CfA 


[H ens cgse] 








[HJ [50], Palomar 










QUEST Survey [EQ 








Stage II 


Pan-STARRSl [62], 


Planck [651 EH 


SDSS II [71], 


Pan- 




HST [S3], KAIT [M] 


M\, SPT [HHJ 


SDSS III 


STARRSl, 






EH], ACT [70] 


[72], BOSS 

LAMOST 
[75], WEAVE 
[76] 


DLS [771 EH], 
KIDS [79] 


Stage III 


DES [HQ], Pan- 


ALPACA, 


DES, HET- 


DES, Pan- 




STARRS4, 


CCAT [H3] 


DEX [Hi], 


STARRS4, 




ALPACA [HI], ODI 




BigBOSS 


ALPACA, 




[82] 




[HH], AL- 
PACA, 
SuMIRe [HS] 


ODI 


Stage IV 


LSST [HZ], WFIRST 


EPIC [HSlEo], 


LSST, 


LSST, SKA, 




[HH] 


LiteBIRD [9H 


SKA [M], 


WFIRST, 






E2], B-Pol [HH] 


WFIRST, 
Euclid [93] 


Euclid 



Table 1: Dark energy projects. Classification is taken from ref. [32]. Note that the DETF 
report was published in 2006, and thus many Stage II projects are now shifted to Stage I. 

We should mention that there exist many other projects planned to survey dark energy. 
Those tremendous efforts may also help us to understand the origin and nature of dark 
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energy. 



2 Cosmological constant 

For explaining the observed accelerated expansion of the Universe, the simplest solution 
was to borrow Einstein's idea of vacuum energy, namely cosmological constant [71 196] . 

Einstein was seeking statistic solutions (a = 0), so he proposed a modification of his 
equation. Einstein's equation with the constant A is given by 

Rfiv - ^g^R - hg^u = 8nGT^. (l) 

With this modification, the Friedmann equations become 

a 2 k 8nG A . , 

-2+~2=— P+o (2) 



and 

The corresponding action is 

S = 



a 2 a 2 3 3 



a 4nG A 
a S 6 



16^G J d ^ x ^iR + 2A ) + ^ 



where g is the determinant of the metric tensor g^ v and R is the Ricci scalar. And Sm 
denotes the matter action. As we see here, the cosmological constant is a constant term 
in the Lagrange density. For the cosmological constant, p\ = —pA = = const., i.e., 

w A = -1. 

The small positive cosmological constant has been supported by a number of observa- 
tions. Indeed the cosmological constant is a perfect fit to the dark energy data, even if we 
cannot explain it. There are two cosmological constant problems [971 [981 EH]- The first 
one is why the vacuum energy is so small or does not gravitate at all, and the second one 
is why it is comparable to the present mass density. 



2.1 Fine-tuning problem 

If the origin of A is a vacuum energy, there is a serious problem of its energy scale. The 
vacuum means a state of minimum energy. For example, for a harmonic oscillating particle 
the potential is of the form V(x) = \uj 2 x 1 . In this case, the vacuum state is given by the 
particle sitting motionless in the minimum of the potential, i.e., x = 0. But quantum 
mechanically we cannot obtain position and momentum of the particle at the same time 
because of the uncertainty principle. Instead, we know that the lowest energy state has 
an energy E = \ftw. 
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We could consider a quantum field as an infinite set of harmonic oscillators, where the 
minimum energy of such a field should be also infinite. However, if we trust our theory 
only up to a certain cutoff, like the Planck scale M p i, from dimensional consideration we 
would get the form 

PA ~ ftMj. (5) 

Indeed, we measure the vacuum fluctuations by the Casimir effect in the laboratory. 
For the Planck scale M pl = (87TG)- 1 / 2 ~ 10 18 GW, we expect 

p A ~ (10 18 GeV0 4 ~ 2 x 10 110 erg/cm 3 . (6) 

However, most of cosmological observations imply 

P°a s < {l^ 12 GeVf ~2x 1(T 10 'erg /cm 3 . (7) 

There is 120 orders of magnitude difference between the theoretical expectation and the 
observational value. This discrepancy has been called 'the worst theoretical prediction in 
the history of physics!'. 

Why is the vacuum energy so small? Does it somehow cancel out exactly a factor of 
10 120 ? These are major outstanding issues in physics and cosmology. Some supersym- 
metric theories predict a cosmological constant that is exactly zero. In supersymmetric 
theories, the number of fermionic and bosonic degrees of freedom are equal. The energy 
of the vacuum fluctuations per degree of freedom is the same in magnitude but opposite 
in sign for fermions and bosons of the same mass. Therefore the fermion and boson con- 
tributions cancel each other and the total vacuum energy density vanishes [SJ. But it is 
not helpful to solve the problem because supersymmetry has to be broken today, as it is 
not observed in nature. 

If supersymmetry is broken, supersymmetric partners can have different masses of order 
^■susyi where Asusy is the supersymmetry breaking scale. Assuming A S usy ~ ITeV, the 
vacuum energy density becomes ~ (10 3 GeV) 4 which is still 60 orders of magnitude 
larger than the observational limit. 

2.2 Coincidence problem 

The second cosmological constant problem is why p\ is not only small but also of the 
same order of magnitude as the present mass density of the Universe. In other words, why 
does cosmic acceleration happen to begin right now and not at some point in the past or 
future? 

If the vacuum energy were big and dominant from the earlier epoch, there would 
be no chance to form structures in the Universe, like galaxies, stars, planets and us, 
intelligent lives. In other words, observers will only observe the states which are allowed 
for observers. This consideration is the so-called anthropic principle. According to the 
anthropic principle, the Universe may not be determined directly by one specific process, 
but there are different expanding regions at different times and space or of different terms 
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in the wave function of the Universe. Therefore, the value of pa that we measure should 
be just suitable for the evolution of intelligent lives. 

This anthropic consideration seems to give us the explanation for both of the two 
cosmological constant problems, why it is small and why acceleration starts now. Actually, 
we sometimes use the anthropic principle as a common sense. Why does the Earth lie 
exactly on the habitability zone, the narrow range of distance from the Sun at which the 
surface temperature allows the existence of liquid water? The answer could be that, if the 
Earth was not there, we would not be here. 

But still, many scientists are not satisfied with this explanation and are disappointed 
by that they could not explain directly the observed Universe from first principles. 

We are living in our Universe, so we cannot verify whether the anthropic principle 
solves the cosmological constant problem or not. The anthropic explanation of the value 
Pa makes sense only if there is a multiverse with a lot of big bangs and different values of 
Pa- 



3 Quintessence 

The Einstein equation G^ u = 9>nGT^ v determines the dynamics of the Universe. (See 



appendix (A.l) for basics.) From the cosmological constant problem we are motivated to 
find an alternative explanation of dark energy. By modifying the left hand side of the 
Einstein's equation, we get the modified gravity models. By modifying the right hand 
side, we get the modified matter models. The idea of modified matter models is that the 
energy momentum tensor T^ u contains an exotic matter, which provides negative pressure. 
The following four sections are devoted to this possibility. 

The first suggestion to solve the cosmological constant problems is quintessence [9"9l 
11001 11011 1102j . The name quintessence means the fifth element, other than baryons, dark 
matter, radiation and spatial curvature. This fifth element is the missing cosmic energy 
density component with negative pressure which we are searching for today. 

The basic idea of quintessence is that dark energy is in the form of a time varying 
scalar field which is slowly rolling down toward its potential minimum. The full action 
including quintessence is given by 



d 4 Xy/—g 



+ S 



167rG T 

The evolution of the scalar field is governed by 

4> + 3H<p + V\4>) = 0, (9) 

where <fi is assumed to be spatially homogeneous. 

The energy density and the pressure of the scalar field are 

pq = \^ + v{4>), pq = ^ 2 -v(0), (io) 
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where ^(j) 2 is the kinetic energy and V{4>) is the potential energy. 
We define the equation of state parameter w as 

PQ |0 2 - V(<j>) 

WQ = = =r~. 11 

Q pq \4> 2 + v^) 

which has a range of — 1 < w < 1. In this range, we are interested in the negative pressure: 
— 1 < w < 0. If the scalar field evolves very slowly so that the kinetic energy term is much 
smaller than the potential energy term, then w is close to —1 and the scalar field behaves 
just like the cosmological constant. 



3.1 Thawing or freezing 

In the equation of motion ([9]), the second term acts as a friction term. By varying the 
friction term, the quintessence model can be dynamically classified [103J. 

• Thawing model 

In the thawing model, the field at early times has been frozen by the Hubble friction 
term 3H(p, and it acts as vacuum energy. When the expansion rate drops below 
H < \/V" {(/)), i.e. underdamped, then the field starts to roll down to the minimum 
and w evolves away from -1. Examples of this model are V(4>) = M A ~ a <j) a for a > 0, 
and V{(j)) = M 4 exp(-/30/M P ) for f3 < ^24^ [TU3] . 

• Freezing model 

If the field is already rolling towards its potential minimum and slowing down, we 
have H > y/V" (</)). In this case, the field is overdamped and approximately constant. 
Examples of this model are V((f>) = M 4+Q </r Q and V((f>) = M i+a (p- a exp(70 2 / M p) 
for a > [103J. In the first case the scalar field energy density tracks that of the 
dominant component (radiation, matter) at early times and eventually dominates at 
late times. Therefore, the tracker behavior partly solves the coincidence problem. 



3.2 Tracker solution 

In the tracker solution |102ll98] . the quintessence component tracks the background density 
for most of the history of the Universe, then only recently grows to dominate the energy 
density and drives the Universe into a period of accelerated expansion. 
The simplest form of the tracker solution is obtained from 

V(4>) = M i+a <f)- a , (12) 
where a > and the value of M is fixed by the measured value of Q m . 
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Early epoch (radiation and matter dominant epochs) 

If the scalar field takes a value much less than the Planck mass in the beginning, its 
energy density is initially given by pq <C pm, where the background energy density 



Pm = Pm + p r - Using the equation of motion for the scalar field with the Eq.(12) 
we have a solution of the field <f>(t) which initially increases as [98] 



2 

4>(t) = 0o ( I ) ^ , (13) 



to 



where we define a dimensionless quantity t = -j-. Then, pq = ^<jr + V(0) decreases 

-2a 

as t 2 + a : 

1 • 2 2a . , 2a —2a 

p Q = -(f) . £2+^ + M A+a t^ oc (14) 
Now the energy density of matter and radiation pm is decreasing faster than the 



Eq.(14) since 



PM = Pm + Pr = Pmfi ( — J + Prfi f ~ J OC t 2 , (15) 

\a m J \a r J 

where the scale factor in the matter dominated Universe behaves as a m (t) oc t 2 / 3 and 
in the radiation dominated Universe as a r (t) oc t 1 ^ 2 . Therefore, the matter dominant 
epoch cannot last forever since pq dominates eventually. 

Does the tracker solution really track the background? To derive the relation 
between the equation of state of the field wq and that of the background Wm, we 

2 

use the known properties, pg oc a~ 3 ^ 1+WQ ' and a oc t z[1+w M) \ 

_2 - 2 (1+™q) 

p Q oc £5+^ oc a~ 3(1+w ^ oc t 1+w m , (16) 
where Q denotes quintessence and M does the background. Comparing exponents 



of t in the Eq.(16), we finally get the equation of state parameter for the tracker 
solution as [102J 

W Q w i i I 

1 2 

which is valid so long as pu S> Pq- This solution is plotted in Figj2] 

For a > 1, we show that wm , i-e., ~ | during the radiation dominated 

epoch and ~ during the matter dominated epoch. In general, wq depends on 



both the effective potential Eq.(12) and %. The effect of the background comes 
from the Eq.([9]), where H carries the information of the background. Another re- 
markable feature is that wq decreases to a negative value (negative pressure that 
we want) after the transition from the radiation dominated epoch to the matter 
dominated epoch regardless of the initial value of Wq. 
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Later epoch (quintessence dominant epoch) 

In later epoch, pg becomes relatively larger and at some time catches up with p M , 
then pq decreases more slowly as 



p Q oc t <h-V (18) 

This is because the field value </>(£) increases as 

0(t)oct*fe. (19) 

The expansion rate H now goes as H oc ^fp~g oc y/V(4>) oc so the scale factor 

a(t) grows almost exponentially 

lna(t)(xt^. (20) 

Thus the quintessence field drives the late time accelerated expansion of the Universe. 
In this solution the transition from pu dominance to pQ dominance is supposed to 
take place near the present time so that both pu and pq are now contributing to 
the cosmic expansion. 

One could also confirm directly that wg ~ — 1 for the later epoch by 

2 



\^ - m _ \ (g) - M^rm 



U2±l/^ ! / n x2 /0 Q \ 2 —4-2« „,~-2a 

Jo . 

where the potential term V becomes more important than the kinetic term as t — > oc. 

As the name tracker solution says, the equation of state wq tracks that of the back- 
ground %. When radiation dominates (wm = |), wq < | and pQ decreases less rapidly 
than p r . When matter dominates (% = 0), wq < and pq decreases less rapidly than 
p m . Finally pQ becomes the dominant component and wq — > —1, as the Universe enters to 
the accelerating phase. See Fig|2j At the transitions from radiation to matter domination 
and from matter to quintessence domination, wq changes almost instantaneously. 

The big advantage of the tracker solution is that, thanks to the existence of a crossover, 
the model does not depend sensitively on the initial conditions for the field. However 
different the initial conditions were, the scalar field keeps rolling down slowly and at the 
end their solutions behave similarly. This is the attractor behavior of the tracker solution. 
Actually, the initial ratio of the tracker field pq to the matter density p m can vary by nearly 
100 orders of magnitude and does not affect the cosmic history |102] . So the quintessence 
model says that the dark energy density is so small because the Universe is old. 

Does it mean that the two cosmological constant problems are solved? The answer is 
no. To make the value of pq at which pg w p M close to the present critical density p c0 , 
we need again fine-tuning of potential energy. 
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Figure 2: A plot of energy density vs. redshift using the equation of state for a tracker 
solution. For an illustration, a = 30 in the Eq.(17) is plotted. The dotted line is for 



radiation, the solid line is for matter, and the thick solid line is for quintessence. 



There is a rough estimate of which value the M should have [98] . If we assume that the 
value of the field (p(t) at this crossover is of the order of the Planck mass, M p i = (8nG)~ 1 / 2 , 



the scalar field evolves very slowly so that V{4>) ~ P<t> ~ Pco — g^g- From the Eq.(12) we 
need 

M 4+a « (8rrG)- a/2 p c0 « {MG)- l - a/2 Hl (22) 

However, the quintessence theory gives no explanation why M should have this value. 

Besides, in order for the field to be slowly rolling today, we require y/V"{(j>) ~ H that 
corresponds to an effective mass of <fi. Thus we must have 

~ H ~ l(r 33 eV, (23) 
which is a very small number in particle physics. 



4 K-essence 

The quintessence model assumes the canonical kinetic energy term ^d^(pd^(f) and the po- 
tential energy term V(<fi) in the action. Modifying this canonical kinetic energy term, the 
non-canonical (non-linear) kinetic energy of the scalar field can drive the negative pres- 
sure without the help of potential terms [104J. The non- linear kinetic energy terms are 
thought to be small and usually ignored because the Hubble expansion damps the kinetic 
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energy density over time. But what happens if there is a dynamical attractor solution 
which forces the non-linear terms to remain non-negligible? This is the main idea of the 
k-essence [IDS} [TPS]. 

The full action including a k-essence term is given by 



d A x\/—g 



1 R + p(^,x) 



+ S M , (24) 



16nG 

where X is the canonical kinetic energy of the field, 

X = ^(V0) 2 , (25) 

and the Lagrangian p(<f),X) plays a role of the pressure px- Here we consider a model 
with 

^=p(0,X)=p(X)/0 2 , (26) 

which has the desired property for dark energy 

For small X, p(X) could be expanded as p(X) = const. + X + 0(X 2 ). If we ignore 
the non-linear term 0(X 2 ) and take an additional potential, then we come back to the 
quintessence model. The scalar field for which these higher order kinetic energy terms 
play an essential role is k-essence. 

The energy density of the k-field is 

p K = (2Xp, x -p)/<p 2 = p/tf (27) 
so that the equation of state parameter for the k-field is 

_Pk P P , OQ , 

p K p 2Xp iX - P 

where t x means derivative with respect to X. If the Lagrangian p satisfies the condition 
Xp t x <C p for some range of X and 4>, then the equation of state is p ~ — p so that we 
have an accelerated expansion solution. We shell find the form of p(X) to satisfy this 
requirment. 

The effective speed of sound c s of k-essence is defined by [107] 

d = M = M. (29) 
P,x P,x 

From observations we know that our Universe is almost flat, so we ignore the curvature 
term in the Friedmann equation : 

H* = tf=*^{p M + pK ), N = lna. (30) 
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-3pi(l + Wi), 



(31) 



The energy conservation equations for the k-essence (i = K) and the background 
(i = M) components are 

dpi 
dN 

where Wi is the equation of state for the corresponding matter (dust and radiation) or k- 
essence. Plugging the energy density (27) in the conservation equation (31 ) and considering 
a homogeneous field 0, we get 



dX 
dN 



P 



P,x 



3(1 +w K ) 



£-1 



/2X 



(32) 



Now we consider the conditions for the stability of the field. We require positive energy 
density, 

p = 2Xp tX -p>0, (33) 
(? s > and the function p(X) to increase monotonically with X, 



p tX = 2Xp : xx+P,x > 0. 



(34) 

It is convenient to use a new variable y = 1/ y/X and to rewrite the pressure and the 
energy density of the k-field as 



pk = g{y)/{4> 2 y), pk = -g'/<J) 2 - 

In this case the equation of state is 

w K = -g I (yg) 

and the sound speed is 

2 p'k 9- g'y 



P' K 



g"y 2 



(35) 

(36) 
(37) 



where prime denotes derivative with respect to y. 

Using this new variable, the stability conditions can be expressed as 



and 



p = -g'>o 
p,x = \y z g" > o. 



(3f 



(39) 



Equally g' < and g" > so that g is a decreasing convex function of y. See the Figure 
4.2. 

Rewriting the equation (32) in terms of the new variable, we get the equation of motion 
for the k-field, 



dy = 3 (w K (y) - 1) 
dN ~ 2 r'(y) 



r(y) - 



Pk 

Ptot 



(40) 
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Figure 3: A plot of g(y). The points indicate the different attractor solutions. R is for the 
attractor solution during the radiation dominated epoch, S is the de Sitter attractor at 
the onset of matter dominated epoch and K is the attractor as k-essence dominates. The 
tangent to the curve at the de Sitter attractor goes through the origin. From Ref.([105j) 



where 



9 , 



1/2 



3 (g-g'y) 



From (31 ) we get 



d(p K Ptot) q Pk A Pk , . . . , 

-n~r = 3 — 1 ){w M - w K (y)). 

dN p tot \ ptot. 



(41) 



(42) 



These equations of motion describe the dynamics of the k-essence model. 



4.1 Attractor solutions 

If every point that is close to some point A is attracted to A, we call the point A an 
attractor. The attractor solutions for k-essence are classified into two types. The first 
one is tracker solution in which k-essence mimics the equation of state of the background 
component in the Universe. In the second case k-essence is attracted to an equation of 
state which is different from matter or radiation. 

We want to have a tracker solution y(N) which satisfies WK{y{N)) = wm = const. 

W K {ytr) = -^\y=y tr = W M (43) 

and Pk/ Ptot is fixed so that the solution satisfies y(N) = y tr = const.. Hence y tr changes 
depending on the epoch and %. The point y tr that satisfies those conditions is the so- 
called attractor solution. Plugging this solution, dy/dN = and d(px/ Ptot)/dN = 0, in 



the equations of motion for k-essence (40), (42) we get a simple relation 



r\y tr ) = ?*- < 1, (44) 

Ptot 
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where we considered px < Ptot- l e., in the range of r(y) > 1 there is no attractor solution. 
Let us consider four possible attractors for k-essence in the different epochs. We con- 



sider positions of attractor solutions, characterize attractor solutions using the Eq.(43) 



and check existence of attractor solutions using the Eq.(44). 
• Radiation tracker 

Radiation has positive pressure and the k-essence pressure px is proportional to g(y). 
Thus the radiation attractor should be located in the region g(y) > 0, equally in the 
range y < yo, where yo is the point which satisfies giyo) = 0. See Figure 4.2. The 
radiation attractor satisfies wr^r) — Wr — g- By substituting this in the equation 



(43), we get 

y R g'(y R ) = -Hvr)- (45) 

During the radiation dominated epoch the k-essence energy density is attracted on 
the radiation attractor and fixed by 

— ) =r 2 (y R ) = -2g'(y R )yl (46) 

Ptot J R 



where the Eq. (41) has been used. 

The most likely range for r 2 (y R ) = i s ~ 10 -2 |106j . so that the cosmic 

acceleration begins at the present epoch. If r 2 (y R ) is much smaller than 10 -2 , the 
k-essence energy density at the dust-radiation equality is so small that it would not 
have overtaken the matter density today. On the other hand, if r 2 {y R ) is much 
greater than 10 _1 , the expansion rate in the early Universe would be so huge that it 
would spoil the predictions of the primordial nucleosynthesis [108J. 

Dust tracker 

The dust attractor corresponds to zero pressure so that WKiyo) = »d = and 
givo) — at the dust attractor y = y^. Thus this attractor should locate at the 
point where g(y) goes through zero. For the dust attractor, we have 

(—) =r 2 ( yD ) = - g -g>(y D )yl. (47) 
\PtotJ D 8 

As we discussed above, the dust attractor exists only if r 2 (yn) < 1. 

De Sitter attractor 

Now we consider the radiation-dust equality phase and at that point the k-essence 

energy density is negligibly small, px <C Pm- Then we have r 2 (ys) = ( — ) — > 

as y(N) — >■ ys in the equation (40). Since r oc (1 + Wk), we finally get wk(vs) ~ — 1 
and 

g(ys) = g'(ys)ys, (48) 
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where the tangent of g(y) at ys passes through the origin. See Figure 4.2. Note that 
we can always find such a point ys for a decreasing convex function g(y). Thus the 
de Sitter attractor (S-attractor) is a generic feature of the k-essence models. 

The assumption px ^ Pm at the dust-radiation equality is also required to satisfy 
the nucleosynthesis constraints |108j . In this scenario, the k-essence approaches to 
the de Sitter attractor shortly after the onset of the matter dominate epoch and 
behaves like a cosmo logical constant, wk — > — 1. The k-essence energy density p% 
becomes very small and freezes at some fixed value. 

Since the matter density pu decreases while px remains constant, the k-essence 
density eventually overtakes the matter density. In this regime, the relation p% <C Pm 
does not hold, the k-field moves to the K-attractor as described below. 

K-attractor 



When pk overtakes pu and dominates, — > 1, then according to the equation (|40|) 



it follows r(yK) = 1 as y(N) ~ yx- This solution describes a power-law expanding 



Universe as follows [104J: substituting r(yx) = 1 in the definition of the r(y) Eq.(41 ), 
we get the relation 

. . 2y/2 1 . . 

1 + w K {y K ) = — = const. (49) 

6 V-9kV 2 k 

and it follows 

a oc t 3 ( 1+ W) = tV-9' K y 2 K/ 2 . (50) 

If it satisfies ~g' K y 2 K j1 > 1, then the solution describes power-law inflation. 

Now let us consider the existence of the K-attractor. We restrict the range of K- 
attractor to yD < Uk < Us so that the k-essence dominant Universe does not have 
positive pressure. If r{y D ) > 1, i.e., there is no dust attractor but exists S-attractor 
V(ys) — 0], then the K-attractor [r(y^) = 1] must exist between y D < y K < y s 
since r{y) is a continuous function. In this interval y > yo, it has negative pressure 
[9(Uk) < 0] and induces the power-law cosmic acceleration. On the other hand if 
t^Ud) < 1, i-e., there is a dust attractor, then since r\y) < for y > yp, where 



3 g"y 



r 



4\/2 v /= ^ 



j(w K -l), (51) 



there is no point y = yx > yD where r(y) = 1. It means, there is no K-attractor at 

vd <y < ys- 

As we considered above, there are two possible scenarios; one is without dust attractor 
(R — y S — > K-attractor) and the other one is with dust attractor (R — > S — > D-attractor). 
According to the first scenario, k-essence is attracted to y = y^ in the radiation dominated 
epoch. At matter dominated epoch, p K drops sharply by several orders of magnitude and 
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k-essence skips the point y = yo- As y ~ ys, Pk freezes and overtakes p m . And then y 
relaxes towards yx- In this scenario, our current Universe lies on the transition from ys 
to y K . 

For the second scenario, with dust attractor [r 2 {%)£,) < 1], we consider again two sit- 



the cosmic acceleration. In the case r 2 (yo) = yff^J ~~ 1> the expansion of the Universe 
would be accelerated before the k-field reaches the dust attractor. In this case, the k- 
essence approaches first the S-attractor, freezes for a finite time, is attracted towards the 
dust attractor, and the Universe decelerates its expansion. This scenario is called late dust 
tracker because the dust attractor is reached long after the matter domination has begun. 
According to the late dust tracker scenario, the cosmic expansion returns to pressureless, 
unaccelerated in the long-term future. 

4.2 Comparing to Quintessence 

Our goal is to have a model which solves the cosmological constant problems. In that 
sense, the dynamical attractor behavior of both the quintessence and the k-essence model 
have a big advantage. The cosmic evolution in this model is insensitive to initial conditions 
because the k-field is attracted to the attractor solution wherever it started. Moreover it 
solves the coincidence problem by explaining why the cosmic acceleration is started at 
such a late stage shortly after the onset of the matter dominated phase. However, both of 
the quintessence and the k-essence do not solve vacuum energy problem. 

Comparing to the tracker solution in the quintessence model, the quintessence field 
tracks the radiation and matter background, and needs a potential energy fine-tuning at 
the quintessence-matter crossover stage. But the k-essence field tracks only the radiation 
background (for no D-attractor scenario), and does not need a potential energy term thus 
it is free from fine-tuning that arose in quintessence. 

For no D-attractor scenario, Wk is increasing today from —1 towards its asymptotic 
value at K-attractor. The numerical value of the effective equation of state in this case is 
w e ff = —0.84 |105j . According to the quintessential tracker solution, the current Universe 
undergoes a phase from w « to w = - 1 and the effective equation of state value 
is w e ff ~ —0.75 |109j . The supernovae observation data are more consistent with the 
k-essence model. 

On the other hand, during the transition from the R-attractor to the S-attractor, there 
is a phase in which wk > 1, i.e., the dominant energy condition px > \pk\ is violated |106j . 
It means, the k-essence energy can travel with superluminal speeds. However, there are 
perfectly Lorentz-invariant theories with the non-standard kinetic term which allows the 
presence of superluminal speeds. Indeed, studies show that in spite of the superluminal 
propagation, the causal paradox does not arise in these theories and in this sense they are 
not less safe than General Relativity |110[ 1111] . 




D 



<C 1, the k-essence could neither dominate today nor cause 
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5 Coupled dark energy and matter 



There is another attempt to solve the coincidence problem which is called coupled dark 
energy [1 121 11131 11141 1115[ 11161 1117] . The dark energy density is of the same order as the 
dark matter energy density in the present Universe. Thus one could imagine that there is 
some connection between dark energy and dark matter. 

The interaction between dark matter and dark energy in form of scalar field is described 
by following modified energy conservation equations 



where 5 is an energy exchange term in the dark sector. The notation m is for dark matter 
and 4> is for dark energy. 

Thus in brief, constructing the coupled dark energy model is finding an appropriate 
form of the coupling 5. There are two major examples [1 131 II 15 j : 



where Q and a are dimensionless constants. We use k 2 = 8nG in this section. Let us call 
those examples as coupling type 1 and coupling type 2. 

5.1 Coupling Type 1 

The first example is a coupled quintessence model assuming an exponential potential and 
linear coupling |113j . The coupled quintessence scalar field equation is 




(52) 
(53) 





5 = KQp m <j>, 



(54) 
(55) 



<P + 3H(f) + 



kQp, 



(56) 



which is equivalent to p^ + 3if(p^ + p^) 



—K,Qp m (j). Here the potential is adopted by 



V{<j>) = V e 



(57) 



The matter component behaves as 





so that matter evolves as 



Pm 



p m0 a e 



-3 k,Q4> 



(59) 



The Friedmann equations for quintessence models are given by 




(60) 
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k 2 ( ■ 4 \ 

H = ~— U 2 + Pm +-Pr) ■ (61) 

We will use the following dimensionless variables |118[ [3J H] 

x=--4- v=-xl^ z = -M (62) 

x 2 , y 2 and z 2 denote the energy density fraction carried by field kinetic energy, field 
potential energy and radiation, respectively. Thus, the density parameters are as follows 

£1$ = x 2 + y 2 , f2 r = z 2 , Vt m = 1 — x 2 — y 2 — z 2 . (63) 



From the equations above (56), (60) and (61) one gets 

x' = -3x + ^\y 2 + \x{3 + 3x 2 - 3y 2 + z 2 ) - ^-Q(l - x 2 - y 2 - z 2 ), (64) 

V = ~\xy + l -y{3 + 3x 2 - 3y 2 + z 2 ), (65) 
z' = -2z + -z(3 + 3x 2 - 3y 2 + z 2 ), (66) 

where ' = -jf- and A = 

a ma kv 

Now we can find the critical points (attractors) that satisfy x' = y' = z' = on which 
the scalar field equation of state is 

x 2 — y 2 

w<h = -r. = const.. (67) 

x l + y l 

The effective equation of state w e ff = Ptot/Ptot is 

z 2 

w eff = x 2 - y 2 + — = VL^w^ + Q r w r (68) 

and w e ff < — | is enough for cosmic acceleration. 

There exist only two attractors which allow accelerated cosmic expansion. However, 
one of them does not contain the matter dominated era and thus fails to explain the large 
scale structure. So we consider just one attractor for this model. According to this scenario 
during the matter dominated era the scalar field has a finite and almost constant energy 
density. This field-matter-dominated era between the radiation era and the accelerated 
era is called (pMDE |113j . The ipMDE is characterized by 

l \ ( ^ n n\ o 2Q2 i 2Q2 teen 

\x,y,z)=\ —,0,0 , i2 = — -, w^ = l, w eff = —-. (69) 
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The (pMDE is responsible for most of the differences with respect to the uncoupled 
quintessence model. For example, the evolution of the scale factor during the <pMDE 
is given by 

aoct^W. (70) 
Finally the field falls into the attractor which is characterized by 

A 2 



{x,y,z) 



A 

71' 



1 



A2 
~6 



.0 



a, 



- h 



A 2 



w, 



eff 



1. 



The scale factor evolves as 



a cx 



(71) 



(72) 



This attractor causes accelerated expansion of the Universe for A 2 < 2. Once this attractor 
is reached, matter density becomes zero. Hence according to this scenario the attractor is 
not yet reached at the present time, but the expansion is already accelerated. 

The CMB data constrain the dimensionless coupling constant to be \Q\ < 0.1 [113J. 



5.2 Coupling Type 2 

In coupling type 1 they assumed a specific potential and coupling from the outset. In 
coupling type 2, interaction potential and coupling structure are determined from the 
requirement — = const. [115] . 



The coupling equation (|53|) is equivalent to 

- + 3H<j) + V 4 



-5, 



where ^ denotes derivative respect to (p. Defining the coupling as 

8 = -3ffll m EE 3HIU, 



(73) 



(74) 



where IT m = —H^ is a relation between effective pressures, the equations (52) and (53) 
become 



p m + 3H(p m + n m ) = 0, 
+ 3H( P(j> + P<p + n ) = 0. 



(75) 
(76) 



Now we consider the requirement for an attractor solution. The time evolution of the 
ratio p m j is 



Pr, 



Pr, 



P<t> 



(77) 



Defining ^ = ^j^- = ^ and p ee p m + we get 



-3R- 



1 - 14, + ^—n m 

PmPtf) 



(7f 
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We can easily find a stationary solution 



n, 



-IX, 



(70 " 1) 



(79) 



which satisfies ( — 

p<t> 



0. 



Substituting the solution Eq.(|79|) in the Eq.(74) yields 

P<j>Pm 



5 = -3H(^-lY 



P 



-3# (t^ - 1; 



r + f 



JO) 



where r = — = const.. 

p<t> 

Now we consider the stability of this stationary solution against small perturbation, 



Pr, 



Pr. 



st 



Then from the Eq.(|78|) we get 

Pm 



e = 3H 



st 



p n r , 



P<f> Pcj> Pm _ 



3H 



st 



Pd> 



+ e 



st 



(81) 



Pr, 



12) 

We choose Il m = — cp where c > so that the interaction is symmetric in p m and p^. Up 
to the first order in e we have 

r 2 -l 

33) 



e = 3Hc- 



so the stationary solution is stable for r < 1. From the Eqs.(74) and (80), the value of the 
constant c is 

c = r 1 - l4> (84) 
(1 + r) 2 1 ] 

and the constant c has positive value for 7^ < 1. Since p today, the stability condition 
which follows from the first expression in ( 82 ) corresponds to 



P Pr, 



< 0. 



55) 



We seek a solution for negative pressure which implies 7^ < 1. Combining this with the 

n 

in, 



Eqs.(74) and (80), we have conditions : H m < 0, S < which lead 

\P\ 



< 



Pr, 



P 



, "><t>+ r 



Connecting the Eqs.(52), (53) with the Eq.(80), we get p m oc p^ oc a~ 3 ~^+ r . Putting it 
in the deceleration parameter (189) implies that the power law accelerated expansion will 
occur for 

r + 3 70 < 2. (87) 
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Combining energy density evolution related to scale factor and the Friedmann equation, 
one gets the time evolution of the scale factor and the field energy density as well. Using 
the definition of 7^, we find finally 

1 7<A 1 + r 1 



6ttG V 2 7 (7^ + r) 2 t 2 



and V((f>) t<j> = -XV ((f>), where A = ^^(7^ + r) so that 

V(4>) = V e- x V-*°\ (89) 

See Ref. |115j for detailed calculation. 

Defining an interaction potential Vmt^ = |, one gets 

V M = ~ 2 1^~ 7 ^ )- (90) 
7^ + r 2 - 70 

Introducing an effective potential 

V e// = V(0) + V^ n4 , (91) 



the coupling equation ( 73 ) becomes 

+ 3#0 + K//,0 = O. (92) 



Considering the condition for accelerated expansion (87), A is restricted to 



A 2 < 247 ^ 7 f )2 . (93) 
(1 + ^)70 

In this model the potential is not an input but derived from the coupling that satisfies 
required property of the attractor solution. But the explanation for how the interaction 
is exactly started at the transition era from decelerated to accelerated expansion is still 
missing. 



5.3 Chameleon mechanism 

As mentioned in section |3j the quintessence field mass must be of order of H . If we 
consider it as a coupled field with matter and assume the interaction is as strong as 
gravity, the coupling must be tuned to a small value (\Q\ < 0.1) to satisfy the test of 
the equivalence principle. The equivalence principle provides that, gravitational mass and 
inertial mass are the same and the laws of gravity are the same in any inertial frame. 

There is a suggestion which allows scalar fields to have couplings to matter of order 
unity. In this scenario the mass of the scalar field depends on the local matter density. 
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In high density environment like on the Earth the field is massive, but in the low density 
environment like in the solar system the field is essentially free. Such a scalar field is 
named chameleon |119[ I12U] as it changes properties to fit its surroundings. 
The action of a chameleon scalar field is given by 



d A x\/—g 



1 -R+lg^d^-V^) 



16nG T 



+ d 4 xC m (^,9Hl), (94) 



where the first integration term is the same as the normal quintessence action and the 
second term is a matter term which is coupled to the chameleon scalar field. Each matter 
field ipm couples to a metric gfy which is related to the Einstein frame metric g^ v by a 
conformal transformation 

9$ = e 2KA V, (95) 

where are dimensionless constants. So the scalar field interacts with matter through 
a conformal coupling in form of e Kl3t ^. 

The field potential is assumed to be of the runaway form and satisfies V — > as — > oo 
and |V^,| — > oo as — > such as 

V(<p) = M 4+ "0-" (96) 

which is familiar from the subsection 13.21 

Varying the action (94) with respect to 0, a similar calculation to the Eq.(225), we get 
the equation of motion for 

d v d»<j> = -Vj - ^e^g^T^, (97) 

i 

where T$ = (2/ ^—g^)S£ m /Sg^ is the stress-energy tensor for the ith form of matter. 
For non-relativistic matter, we have g^T^} p i} where pi is energy density. Introducing 
the energy density pi = pie 3Kl3i ^ which is conserved in the Einstein frame, the Eq.(97) 
reduces to 

d v d v <j> = -V ><t> - "Pine"**- (98) 

i 

If we define the right hand side of this equation as —V e ff^, the effective potential is dubbed 

F e// (0)=V(0)+J> e ^. (99) 



as 



For monotonically decreasing V((p) and positive the effective potential V e ff has a min- 
imum. See Figure 4.3. 

Calculating the point that satisfies V e ff^ = and K//,</></> = m2 ; we can define m j n as 
the value of at the minimum and its mass m min . Furthermore we recognize here that 
larger values of pi corresponds to smaller min and larger m min , i.e., the denser environment 
follows the more massive chameleon. 
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Figure 4: The chameleon effective potential V e ff (solid curve) is the sum of two contribu- 
tions. One from the actual potential (dashed curve), and the other from its coupling 
to the matter density p (dotted curve). From Ref.( |119] ) 



Let us consider a solution for a compact object which is static, spherically symmetric 
with radius R c , homogeneous density p c and total mass M c = 4iiRlp c /3. Then the Eq. 



(98) becomes 

where all /3, are assumed to be of the same value (5. And p(r) = p c for r < R c , p(r) = p^, 
for r > R c , where p^ denotes the surrounding homogeneous matter density. 

We define C and 0oo as the value of that minimizes V e ff with p = p c and p^, re- 
spectively. And the corresponding masses of small fluctuations are m c and m^o. Boundary 
conditions are given by d<p/dr = at r = so that the solution is non-singular at the 
origin and — > 0oo as r — > oo so that the 0-force on a test particle vanishes at infinity. 

Now we consider a large object. Inside the object the field minimizes V e ff at ~ C . 
Outside the object the field evolves as <fi ~ exp(— m^r^jr and approaches to 0^ for 
r 3> R c . During the transition, at r = R c , and d(j)/dr should be continuous thus at this 
point a thin shell of thickness AR C below the surface is introduced where the field grows. 



Solving the Eq.(lOO) for outside the object, one gets 



0(r) - -0OO (l - [9 ^ C> ) ^— r + 0oo, (101) 

where <3> c = GM C /R C is the Newtonian potential of the object. 

In case of the thin shell condition, AR C /R C 1, and the inside of the object being 



much denser than the outside, C -C 0oo, the exterior solution (|101|) becomes 

MJ \ ^ 
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where AR C /R C can be adjusted as 



AR C _ /tQoo - <f> c ) 

i? c ~ 6/3$ c • (103) 

In case of small objects, in the sense of AR C /R C > 1, however, the thin shell condition 
is not satisfied. Since their entire volume contributes to the <f)-fie\d outside, the fraction 
of the infinitesimal volume in the thin shell to the total volume 3AR C /R C becomes 1, so 
the exterior solution is 

m - - [-£) ^— r — + (104) 

Consider the potential V(<p) = M 4+n cj)~ n , take n — (3 — 1 and apply these to the 
Earth which must have a thin shell. It follows that the interaction range (which is anti- 
proportional to the field mass) is of order 1mm on the Earth and of order 10 — 10 4 AU in 
the solar system |119j . 

The chameleon with thin shell effect satisfies tests of gravity in laboratory as well as 
from the solar system data |119j . Moreover the chameleon mechanism predicts that the 
magnitude of equivalence principle violations and fifth force are much greater in space 
than on the Earth. The chameleon force on a test particle of mass M and coupling (3 is 
given by 

= -K/3MV0, (105) 
so the (f) can be thought as a potential for the fifth force. 



get the fifth force on a test mass M and coupling of magnitude 

, : K 2 (3f3 t f:\AR \ M©M 



On the Earth, calculating V0 from the Eq. (]102|) and substituting in the Eq.(105) we 

4, V « 9 I ' (1 ° 6) 

with the Earth radius R® and the Earth mass M®. The Earth has a thin shell , Ai? e /i? e 
1, hence the fifth force on the Earth is suppressed. 



But in satellite in space, calculating the V0 from the Eq.(104) without thin shell and 



substituting in the Eq.(|105|), we get the significant fifth force between two bodies of mass 

, * (107) 



Mi and M 2 and its coupling (3 t and (3 2 



Then the total force, i.e., the gravitational force plus the chameleon- mediated fifth force 
between two masses is given by |119j 

fl_0^ (1 + JAA) , (10 8) 
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hence, the effective Newton's constant is G e // = G(l + 2/3i/3 2 ) • 

In a similar way chameleon predicts an equivalence principle violat ion. We can read 
out the extra acceleration component a<p = — from the Eq.(105). The Eotvos pa- 
rameter which denotes relative difference in free-fall acceleration for two bodies of different 
composition is defined by 

V^ — = ^, (109) 

where a at is the Newtonian acceleration on the Earth. If future gravity tests in satellite 
measure an effective Newton's constant which differs by order unity from the measured 
value on the Earth, or find an equivalence principle violating signal stronger than allowed 
by a laboratory experiment, as the chameleon cosmology predicts, it will be a strong 
candidate for dark energy. 



6 Unified dark energy and matter 



6.1 Chaplygin gas 

To avoid the potential energy fine-tuning of the quintessence, instead of the form of po- 
tential, we take a change in the equation of state of the background fluid |121] 1122] . This 
exotic background fluid, the so-called Chaplygin gas, has the following equation of state: 



p = — 



A 



(110) 



where A is a positive constant. Inserting this equation of state into the energy conservation 



equation (187), we obtain 



p{t) 



A 



B 



3(l+a) 



a 



i 

l+a 



(111) 



where B is an integration constant. 

In the early epoch, a C 1, the Chaplygin gas energy density behaves as p oc a -3 which 
corresponds to the matter dominated Universe. In the late epoch, a > 1, the energy 
density behaves as p ~ A l ^ a+a ^> = const, which corresponds to the de Sitter Universe. 
Thus the single fluid, Chaplygin gas, behaves as dark matter in the early epoch and dark 
energy in the later epoch. That is why we call it a unified model of dark energy and dark 
matter, and it is thought to be an attractive feature that explains both dark sectors in 
terms of a single component. 

The effective sound speed for the Chaplygin gas is given by 



dp 
dp 



—aw, 



(112) 
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where w is the equation of state parameter of the Chaplygin gas. Substituting above 
Eqs.flTTo| and (fill]) into the Eq.([TT2|)we obtain 



a: 



1 + 



B/A ' 



(113) 



In the Chaplygin gas model the sound speed is small at the early epoch and becomes 
larger at the late epoch. The problem is, the large sound speed at the late epoch leads 
to growth of inhomogeneities. If so, we should have observed fluctuations or blow up in 
the matter power spectrum in the large-scale structures, but they are not observed. The 
observation constricts the upper bound on the values of a |123] 



< 10" 



(114) 



Thus, the Chaplygin gas model with \a\ ^> 10 -5 is ruled out. Moreover the Chaplygin gas 
mode is indistinguishable from the ACDM model which implies a = — 1 and w = const. = 
in the early epoch and —1 in the later epoch, and c s = —aw = 1 at the present epoch. 



6.2 K-essence as unified dark energy 

The problem of the Chaplygin gas model can be avoided by constructing unified model of 
dark energy and dark matter using the k-essence model |124j . 

For the Lagrangian density of k-essence, p(X) as in the Eq.(24), we assume that p(X) 
can be expanded about its extremum at X = Xq in the form 



p(X)=p +p 2 {X-X 



(115) 



where po and p 2 are constants. The k-essence pressure Pk = P and the energy density 
Pk = 2Xp 5 x —p satisfy the energy conservation equation pk + 3H(pk +Pk) = and follow 



(p,x + 2Xp :XX )X + 6Hp jX X = 0. 
Around the solution X = Xq we introduce small perturbation e = 



(X-X ) 

x 



(116) 
< 1. Sub- 



stituting the Eq.(115) into the Eq.(116) and calculating up to linear order we obtain 

e = -3He. 



From e = X / ' Xq and the Eq.(117) we obtain a solution 

X = X [l + e 1 (a/a 1 )- 3 ] 



(117) 



(118) 



in terms of new constants e\ and a>\. X w Xq when ei(a/ai)~ 3 1. 

Plugging the Eq.(115) into px = 2Xpx —p and using the solution (118) around 
X = X , we get 

Pk ?t po, Pk ^ ~Po + 4p 2 ^o e i( a / a i)~ 3 - ( 119 ) 
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Now we obtain the equation of state for k-essence as unified dark matter 

-l 

4p 2 .,9 ( a\ 

w K ~ 



1 - -Z±Xfa 



(120) 

It follows wk — > at the early epoch and wk — )■ — 1 at the late epoch. 



The resulting effective speed of sound as defined in the Eq.(29) is 
thus the large sound speed problem is avoided. 



<^hJ a ) «i- (i2i; 



7 f(R) gravity 



For modified matter models we have considered a modification of the right hand side of 
the Einstein's equation G^ u = 8ttGT^ u . Now we consider to modify the left hand side of 
the Einstein's equation. It means we do not search for matter causing cosmic acceleration 
anymore, but rather modify gravity itself. Therefore in the action we do not add a scalar 



field term as we have done in the Eq.([8j) and the Eq.(24), but rather modify the Einstein- 
Hilbert action term. 

To modify the Einstein-Hilbert action one may consider an action with higher-orders 
of curvature invariants, for example, 

S = --^ [ d A x^(R + aR 2 + pRuuR^ + -?R 3 + ■■■)■ (122) 
1o7tG J 

The higher-order terms can be either of fundamental origin or they can arise as a result 
of vacuum polarization [9]. In general Einstein gravity causes second order equations 
of motion. But every modification of the Einstein gravity introduces higher-derivative 
terms. I.e, in addition to the gravitational waves, the gravitational field has extra degrees 
of freedom including a scalar field [9j. 

7.1 f(R) cosmology 

Consider an action 

1 



167rG 



d 4 x^f(R) + / fx^M, (123) 



where f(R) is an arbitrary function of the scalar curvature R. Using a conformal transfor- 
mation g^y —7- g^y = (df '/ l dR)g flu , it is possible to show that the higher derivative gravity 
theory (Jordan frame) is conformally equivalent to Einstein gravity with an extra scalar 
field (Einstein frame) \l'25\ [pj. We call this modified gravity theory as f(R) gravity 
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[125[ I126[ I127j . If the scalar field potential satisfies the slow-roll condition, then we have 
an inflationary solution. Thus f(R) gravity causes inflation due to the potential of a scalar 
field. 

Varying the action (123) with respect to the metric g^, it yields the modified Einstein 
equation 



-Jg^u - Vf.Vuf^R + g^ n f,R = 8-trGT^ 



(124) 



where jR denotes the derivative with respect to R, V M is the covariant derivative and □ is 
the D'Alembert operator defined as □ = g^ v V ^V v . T^ u is defined as in the Eq.(207). The 



trace of the Eq.(124) is given by 

3nf >R + f ;R R-2f = 87rGT, 



(125) 



where T = g^T^ = p M - 3p M - 

To get an inflationary solution we consider the de Sitter point. The de Sitter space is 
a vacuum solution (T = 0) with constant positive curvature, i.e., R = const.. Thus at this 
point Of R = so that 

f, R R-2f = 0. (126) 

If a dark energy model based on the f(R) gravity satisfies this condition, the late-time de 
Sitter solution can be realized. 

Assuming a flat FLRW metric, we get modified Friedmann equations from the Eqs. ( 124 ) 
and ([125} [125] : 

3f, R H 2 = 8nG Pm + (f >R R - f)/2 - 3^, (127) 
2f >R H = -8nG Pm - f, R + Hf iR , (128) 

where the dot means derivative with respect to the cosmic time t and the radiation com- 
ponent is ignored. To compare f(R) gravity with observation we rewrite the Eqs. (127) 
and (128) as 

3 AH 2 = 87rG(p m + p DE ), 
-2 AH = 8nG(p m + p DE + Pde), 

where A is some constant and p m = 0. pde and pde are defined by 

8nGp DE = (f, R R - f)/2 - 3Hf tR + 3H 2 (A - /,*), 
8nGp DE = f >R + 2Hf jR - (f >R R - f)/2 - (3H 2 + 2H)(A - f >R ) 



(129) 
(130) 

(131) 
(132) 



so that the continuity equation (187) holds for pde and pde- From Eqs. (129) and (130) 
we obtain 

p DE 2AH + 3AH 2 

p DE 3AH z - 8nGp m 
In order to recover the standard matter era in the past we can choose A = 1 in the 
Eq.(129). The equation of state Wde can be smaller than -1, i.e., a phantom equation of 
state, before reaching the de Sitter attractor 
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7.2 Cosmological and local gravity constraints 

The modification of gravity affects the large-scale structure as well as local gravity. Study- 
ing them, we can observationally distinguish between the f(R) gravity model and the 
ACDM model. At the same time we obtain constraints that the viable f(R) has to satisfy. 

First of all, we need f\n{R) < to avoid anti-gravity. And the effective scalar field 
mass of the f(R) gravity model is given by |129] 



M f(R) 



3/, 



RR 



(134) 



in the regime M 2 ^ 3> \R\. Thus we need the condition f t R R (R) > to avoid a tachyonic 
instability related with negative mass squared |129j . From those two stability conditions 



and the existence of a late-time de Sitter point given by the Eq.(126), it is stable for 
< ^« < l [T28]. 

The first suggestion is f(R) = —R + aR 2 |13Uj . But this model is not suitable to 
explain current accelerated expansion of the Universe because R 2 <C \R\ at the present 
epoch |128j . The next propose f(R) = -R - a/R n (a > 0,n > 0) |126] is also ruled out 
by the above constraints. 

The models that fulfill all such constraints are for example 



f(R) 
f(R) 
f(R) 



-R — [iR c 



-R — [iR c 



(R/R c 



x 2ll 



(R/R c ) 2n + 1 



—R — fiR c tanh 



1 + 



\R\ 

Rr 



Rl 



with n, fi, R c > 0, 

with n, /i, R c > 0, 



with /x, R c > 0. 



(135) 
(136) 
(137) 



These models satisfy the relation f(R = 0) = 0, i.e., the cosmological constant disappears 
in a flat space time. 

The ACDM model and the f(R) gravity model predict different structure formation 
histories. For the large-scale structure observations we consider the wavenumber k in sub- 
horizon scales k/a 3> H. By quasi-static approximation one gets the equation for matter 
density perturbation 5 m |131j 



5 m + 2H5 m - 4irG eff p m 5 m ~ 0, 
where the effective gravitational constant G e // is defined by [131 J 

G 1 - Amk 2 /(a 2 R) __ Rf >RR 



G 



m 



(138) 



(139) 



f >R l-3mk 2 /{a 2 RY ~ f, R 

Here m characterizes the deviation from the ACDM model which denotes f{R) = —R—2A. 
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If the deviation from ACDM is small, i.e., —mk 2 /(a 2 R) <C 1, then the effective grav- 
itational coupling G e ff is very close to the gravitational constant G so that S m oc t 2 ^ 3 
during the matter dominant epoch. But in the regime —mk 2 /(a 2 R) 3> 1, the effective 
gravitational coupling approaches G e // — — AG/3f t R so that 5 m oc ^(v^- 1 )/ 6 |128j . Later 
we will see that the DGP model implies the same aspect, i.e., on small scales we recover 
the Einstein gravity meanwhile on large scales the modified gravity becomes important. 

Computing the matter power spectrum Ps m = \S m \ 2 , we can justify whether the f(R) 
gravity realizes in nature. Moreover the modified evolution of matter perturbation directly 
affects the shear power spectrum in weak lensing |132j . To be consistent with the local 
gravity constraints in the solar system, the function f(R) needs to be very close to that in 
the ACDM model in high density regions |133j . The SN la data constraint the deviation 
parameter m to be m(z = 0) < 0.3 |128j . In high density regions, \R\ S> \Rq\, the linear 
expansion of R with respect to the cosmological value Rq is no longer valid. In order to 
evade the solar system tests, one derived the non-linear chameleon mechanism |133j . where 
the chameleon mechanism may suppress the fifth forces in such a nonlinear regime. 



8 DGP model 

The brane world model of Dvali, Gabadadze and Porrati (DGP model) |134j suggests that 
cosmic acceleration is a signal of lack of understanding of gravitational interactions. Let us 
suppose that our 4-dimensional (4D) world is a brane which is embedded in 5-dimensional 
(5D) bulk Minkowski space-time with infinitely large extra dimensions. All particles and 
standard model forces are pinned onto the brane world, like dust particles on soap bubbles, 
while gravity is allowed to explore into the 5D bulk. 

8.1 Brane cosmology 

The action for the DGP model is given by 

ML r . M 2 



S= ^ / (FXy/^jjR--^ I d 4 x^g~R + I d*x^g~C m , (140) 



where gAB is the metric in the bulk and R is its Ricci scalar. And is the induced metric 
on the brane and R is the corresponding Ricci scalar. The first and second terms form 
Einstein-Hilbert actions in the 5D bulk and on the brane, respectively. The third term is 
the matter action where C m is the matter Lagrangian confined to the brane. The second 
term may be generated by quantum corrections from the 5D gravity, or its coupling with 
a certain 5D scalar field |134j . 

Capital letters, superscripts and subscripts will be used for 5D quantities (A, B = 
0,1,2,3,5) whereas the 4D coordinates of the brane are x M (/x = 0,1,2,3). The extra 
coordinate will be denoted by y. The induced metric on the brane is given by 

g, u = d,X A d u X B ~g AB , (141) 
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9nu{x) = gn„(x, y = 0). (142) 

Planck mas 

scale r c is denned by 



We define the 4D Planck mass as M p i and the 5D Planck mass as Mtg\. The cross-over 



If the characteristic length scale r = ^x\ + x\ + x\ is much smaller than the cross-over 
scale, r c , gravity behaves as usual 4D theory. In contrast, at large distance gravity slips 
into the bulk so that it occurs a weakening of gravity on the brane. Thus the higher 
dimension plays an important role. 

It can be easily understood through considering the brane as a metal sheet immersed 
in air |135] . Imagine that sound waves represent gravity. If one strikes the metal sheet the 
sound waves propagate along the metal sheet as well as into the air. However the energy of 
the sound wave on the sheet is much denser than in the air. Therefore an observer on the 
sheet does not feel the extra dimension (air) at all. Only after the wave has propagated a 
long distance so that an amount of sound wave energy has been lost into some unknown 
region, the observer becomes aware of the existence of a extra dimension. 

Across the cross-over scale r c , the weak-field gravitational potential behaves as |134l 

V{r)~{ r ~\ l° r r<<rc ' (144) 
yr tor r » r c . 

How does the extra dimension change cosmology on the 4D brane? Let us consider the 
5D space-time metric 

ds 2 = g AB dx A dx B . (145) 
Since we are interested in cosmology on the brane, the 5D line element is given by 

ds 2 = N 2 (t, y)dt 2 - A 2 (t, y)-f ij dx i dx j - B 2 (t, y)dy 2 , (146) 

where 7^ is a maximally symmetric 3D metric and the brane is a hyper-surface defined 
by y = 0. The metric coefficients read |137] 

N(t, y) = 1 + e\y\d(a 2 + k)~ 1/2 , (147) 
A(t,y) = a + e\y\(a 2 + k) 1 / 2 , (148) 
B(t,y) = l, (149) 

where a(t) is the 4D scale factor and e = ±1. Taking y — 0, we obtain the usual 4D FLRW 
form. 

The Einstein equation in 5D bulk is given by [137J 

Gab = Rab - \rqab = 0, (150) 
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where Gab is 5D Einstein tensor. The induced 4D Einstein equation is as follows [137l ll38j 

1 



G 



fin 



-(K llu -Kg flu )=87tGT > 



(151) 



where is the extrinsic curvature on the brane. On the FLRW brane we get the modified 
Friedmann equation \137\ 1136] 



H 2 + 



k 



P 



3M* 



1 e 

4r 2 2r r 



(152) 



where p denotes the total cosmic fluid energy density on the brane which satisfies the 
conservation equation (187). If ^> i.e., in the early Universe, this Friedmann 



equation takes the standard cosmology form: if 2 + \ 



-P- 



For a flat geometry (k = 0) the Eq.( 152) takes the following form 



H 1 



-H 



(153) 



^Y-p- On the other hand, in case of if 1 3> r c the 



In case of if -1 r c the second term of the Eq.(153) is negligible so that we again recover 
the usual Friedmann equation: H 2 
second term becomes important. 

Depending on the sign of e, the cosmological solution Eq.( 153 ) has two different regimes. 
If e = — 1 and H~ l ^> r c , the Friedmann equation approaches 



P 



(154) 



If e = +1, however, the matter dominated Universe (p oc a 3 ) approaches the de Sitter 
solution, 

H^H OQ = -, (155) 

r c 

which causes late-time acceleration. This cosmological solution is very interesting because 
it drives our Universe into self -inflationary regime without dark energy. 

Requiring that the late-time acceleration occurs around the present epoch, we get |136j 



which corresponds to 



Mi 



(5) 



10 - lOOMeV. 



(156) 
(157) 
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8.2 Observational test 

The most important task for observational tests of modified gravity is to distinguish it 
from the ACDM model. In order to discuss about cosmological tests it is convenient to 



rewrite the modified Friedmann equation ( 152 ) in term of the redshift as |136] 



HHz) 



Hl\n k {\ + z f + (v^T + y/n rc + a*(i + ^) 3 ) 2 J, 



(15* 



where the redshift is defined by 1 + z = — . And density parameters are defined by 



__ 8%G p mfi _ -k 



3 flgag' ^"4^' 
We compare this equation with the conventional Friedmann equation 

f 2 (z) = ^ 2 {fi fc (i + zf + n m (i + ^) 3 + n DE (i + 



Setting 2 = in the Eq.(158), we obtain the normalization condition 
which differs from the conventional relation Q k + fi m + Qde — 1- 



(159) 



(160) 



(161) 



For a flat Universe (fi*. = 0) the Eq.(161) becomes 

2 



1 - tt r 



Q r < 1. 



(162) 



Comparing the Eq.(158) and the Eq.(160), we see that f2 rc plays a role of but not 



Different Friedmann equations result in a different luminosity distance d E 



same, as shown in the Eq.(161). And for a flat Universe Vt Tc is always smaller than VLde- 

(163) 



driz] 



l + z 



Sk Hn\/£tt 



dx 



H(x) 



where 



S k (x) 



sin(x) for k — 1, 
x for k — 0, 

sinh(a;) for k = —1. 



(164) 



The apparent magnitude in units of Mpc is given by 

m = M + 5 log di + 25, 



(165) 
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where M is the absolute magnitude. Considering supernovae la as standard candles, M 
is the same for all supernovae, and measuring apparent magnitudes of supernovae we can 
directly compare our model with SN la observations. To use this method we have to know 
the exact value of H . 

In spite of attractive features, the DGP model suffers from observational disfavor and 
ghost instabilities. SN la, BAO and CMB data shows that the modified Friedmann equa- 
tion is less consistent with observations than the usual Friedmann equation, i.e., ACDM 
model |139| . There is a modified version of the DGP model which is characterized by |14Uj 

h-2 H ° 8nG hurt 

H ~ = ~^rP ( 166 ) 



' c 



so that the effective equation of state w e ff evolves from — 1 + -y during the radiation 
dominated era to — 1 + |, during the matter dominated era to —1 in the distant future. 
Recent observational data have constrained a = 0.254 ± 0.153 [139J, thus the flat DGP 
model (a = 1) is ruled out. 

According to the studies of the linear theory about a flat multidimensional space-time 
and a flat brane, the DGP model has a scalar ghost field localized near the brane [14T| . 
However it is possible to get a ghost-free DGP model by embedding our visible 3D brane 
within a 4D brane in a flat 6D bulk |142] . Moreover Galilean gravity also gives rise to the 
possibility of avoiding the ghost problem |143j . 



9 Inhomogeneous LTB model 

The apparent accelerated expansion of the Universe may not caused by dark energy but 
rather by inhomogeneities in the distribution of matter. The basic idea is that, there are 
inhomogeneities on a larger scale, in the form of underdense bubbles, i.e., we live in an 
underdense region of the Universe and describe its behavior, a faster expansion compared 
to the outside, as an apparent cosmic acceleration. 

How do we understand the concept of apparent cosmic acceleration! In a homogeneous 
Universe the expansion rate is a function of time only. But in an inhomogeneous Universe 
the expansion rate varies with both time and space. Thus, if one observes faster expansion 
rate for low redshifts than higher redshifts, in the homogeneous case it is described as 
cosmic acceleration, whereas in the inhomogeneous case it is a result of spatial variation 
with an expansion rate being faster as being closer to us [144] . 

The first underdense void model is proposed by Tomita |145j in form of a local homo- 
geneous void separated from the homogeneous FLRW outside with a singular mass shell. 
There occurs a discontinuous jump at the location of the mass shell. This model is ex- 
tended to a more realistic model with continuous transition between inside and outside of 
the void |144j . In this extended model the Universe is dust dominant and isotropic but in- 
homogeneous. The inhomogeneity is spherically symmetric. This model can be described 
by the Lemaitre-Tolman-Bondi (LTB) spherically symmetric models |1461 11471 f!48j. 
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We have used the FLRW metric under the assumption that our Universe is isotropic and 
homogeneous. Now we renounce the homogeneity, then the line element for a spherically 
symmetric inhomogeneous Universe is given by 



d s 2 = dt 2 - X 2 (r, t)dr 2 - R 2 (r, t)dn 2 . 



(167) 



For the Einstein equation we assume T^ u 
Solving the equation G i = gives 

X(r,t) = 



diag(p, 0, 0, 0), i.e., containing matter only. 
R'(r,t) 



(168) 



yfi+ffr) 

where prime denotes derivative with respect to r and /3(r) is a function of r. The metric 
(167) recovers the FLRW metric by choosing R = a(t)r and (5 = —kr 2 . 

The Einstein equations for the dust dominated LTB Universe are given by |144] 



Hl + 2H\\H 



'R 2 



R 2 
2H\\H\ 



JL 

RR' 

JL 

RR' 



8nGp, 
-8nGp, 



(169) 
(170) 



where dot denotes derivative with respect to t. The transverse Hubble function is defined 
by H± = R/R and the radial Hubble function is defined by H» = R'/R'. Adding Eqs. 
jl69| ) and ( fl70| ), we obtain 

2RR + R 2 = (3, (171) 



and integrating this equation leads to 

H 2 



a (5 
R3 + R 2 ' 



(172) 



where a is a function of r. Thus the dynamical effect of a and (5 are similar to those of 
dust and curvature, respectively. 

We define the deceleration parameter q± = —RR/R 2 . See the Eq.(189). Substituting 
Eqs. (171) and (172) into this expression yields 

Q± 



a 



2a + (3R 



(173) 



Since a > 0, the deceleration parameter has a positive value. Indeed for a spatially flat 
and dust dominated Universe, q± = 0.5. Thus a dust dominated inhomogeneous Universe 
can not be accelerating. 

We define t = as the time when photons decoupled from matter (z ~ 1090) and 
R(r,t = 0) = 0. And we introduce a conformal time rj by dr\ = {^J~f3/R)dt. Integrating 
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Eqs. (169) and (170) with (3 > 0, we get a solution of the Einstein equation for the LTB 



model |144j : 



R = ^(coshr/ - 1) = T^^( cosh V - 1> 



t 



(sinh rj — rf) = m '° , (sinh n — Tj), 



2(3 3 / 2 



(174) 
(175) 



where we chose a = Hj_ Q m fir 3 and (3 = Hj_ Q k>0 r 2 . 

The structure of the underdence void causing apparent acceleration can be expressed 
as [US] 



^m,o(^) ^out iS^in ^out) 

H±fl(r) = H out + (H in — H out ] 



1 — tanh[(r — ro)/2Ar] 

1 + tanh[r /2Ar] 
1 — tanh[(r — ro)/2Ar 
1 + tanh[r /2Ar] 



(176) 
(177) 



where in and out represent quantities inside and outside the void, respectively. Further r 
characterizes the size of the void and Ar is the thickness of the transition shell. 

From the CMB acoustic peak the value of the local density parameter Qi n is constraint 
to be in the range 0.1-0.3, whereas VL out = 1 for asymptotic flatness. The HST obser- 
vations imply H in w 70 km/sec/Mpc, whereas outside the void one requires H out w 50 
km/sec/Mpc to be consistent with VL out . The SN la data constraint r and Ar to be 
r = 2.3 ± 0.9 Gpc and Ar/r > 0.2, respectively [H9] . 



We define relative matter and curvature densities from the Eq.(169) as 



8irGp 



H 2 + 2H»H 



n, 



(178) 
(179) 



The usual definition (197) for the homogeneous case can be recoverd with H± = Hu 



To compare the inhomogeneous LTB model with SN la observations, we need to find 
the luminosity distance. Photons which travel along radial null geodesies, ds 2 = dQ 2 = 0, 
arriving at r = today t follow a path t = t(r) given by 



dt 
dr 



R'(r,t) 



with t(0) = to- The redshift z 



z(r) of photons obeys the differential equation |150] 



(180) 



dz 
dr 



R'(r,t) 



;i8i) 
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with z(0) = 0. Now we have the functions t(z) and r(z). From that we can obtain the 
luminosity distance di given by 

d L (z) = (l + z) 2 R(r,t), (182) 

which is related to the angular diameter distance (Ia(z) = R(r,t). 

The inhomogeneous LTB model matches to the supernovae data and the location of 
the first acoustic peak of CMB temperature power spectrum |144j but it is still challenging 
to reproduce the entire CMB angular power spectrum. The observed isotropy of the CMB 
radiation implies that we must live close to the center of the inhomogeneity. The maximum 
distance to the center is constrained to be smaller than 15 Mpc |151j . I.e., we are located 
at special point in the Universe. This reminds us again of coincidence and the anthropic 
like question, why are we special? 

10 Summary 

We summarize and compare the various dark energy models with respect to attractive 
features, rising problems and possible solutions. 

Cosmological constant 

The simplest explanation for dark energy is the energy associated with the vacuum which 
has constant equation of state parameter w — — 1. The cosmological constant has perfect 
fit to observations as shown in Figure 1.1. But it suffers two cosmological problems. The 
observed dark energy density value is 120 orders of magnitude smaller than the theoretical 
expectation (fine-tuning problem). Moreover dark energy density and dark matter energy 
density are of the same order today (coincidence problem). 

Quintessence 

Dark energy could exist in the form of a scalar field which has time varying equation of 
state. According to the tracker solution, the quintessence component tracks the equation 
of state of the background (radiation in the radiation dominant epoch and matter in the 
matter dominant epoch) and only recently grows to dominate the energy density. The 
tracker behavior allows the quintessence model to be insensitive to initial conditions. But 
it needs fine-tuning of the potential energy as y/V"{(j>) ~ H ~ 10~ 33 eV. Moreover, the 
quintessence model does not solve the first cosmological constant problem. 

K-essence 

In the k-essence model, the non-canonical kinetic energy term in the Lagrangian can 
drive late-time cosmic acceleration without the help of potential energy. The attractor 
behavior allows the k-essence model to be insensitive to initial conditions as well. Unlike 
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the quintessence model, the k-essence field tracks only the radiation background so that it 
is free from fine-tuning which arose in the quintessence model. Moreover the coincidence 
problem is solved by presence of an S-attractor which attracts shortly after the onset of 
the matter dominated phase. However the k-essence model does not explain the smallness 
of the vacuum energy. There is a superluminal sound speed phase, but a causal paradox 
does not arise. 

Coupled dark energy 

Dark energy density and dark matter energy density are of the same order today (coin- 
cidence problem), thus it is natural to imagine that dark energy and dark matter have 
some relation. Coupled dark energy models connect dark matter and dark energy. Hence 
dark sectors exchange energy with each other. But observations constraint the coupling to 
be small (coupling constant \Q\ < 0.1). This constraint may be avoided by adopting the 
chameleon mechanism in which a scalar field has couplings to matter of order unity and 
changes its mass to fit to the local matter density. It is expected that the measurement 
of effective gravitational constant in satellite project test this model. 

Unified models of dark energy and dark matter 

These models explain both dark sectors in terms of a single component by using a Chap- 
lygin gas with exotic equation of state. The Chaplygin gas behaves as dark matter in the 
early epoch and as dark energy in the later epoch. But it has large sound speed at late 
epoch which leads to growth of inhomogeneities so that it is ruled out by observation. The 
large sound speed problem can be avoided by introducing k-essence as unified dark matter 
which provides also w — > at the early epoch and w — > — 1 at the late epoch. 

f(R) gravity 

f(R) gravity modifies Einstein gravity. The higher order gravity terms (in form of function 
f(R)) in the Lagrangian are equivalent to Einstein gravity with an extra scalar field. If 
this scalar field satisfies some condition it can drive cosmic acceleration as well. So the 
f(R) gravity is thought to be equivalent to the quintessence. There are several models 
of function f (R) which fulfill stability conditions as discussed in section [7j Observing the 
matter power spectrum or weak lensing, we can justify whether f(R) gravity realizes in 
nature. The f(R) gravity model is strongly constraint by solar system tests. Adopting the 
non-linear chameleon mechanism, it is possible to build a viable model. 

DGP model 

According to the DGP model, we are living on 4D brane world in a 5D extra dimensional 
bulk. All particles and standard model forces are pinned onto the brane world, only 
gravity can explore into the 5D bulk. It appears in large scale as weakening of gravity on 
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the brane. The modified Friedmann equation in the DGP model allows self-acceleration of 
the Universe. It may explain the weakness of gravity compared to other forces and gives 
a connection to string theory. Moreover the DGP model can solve the first cosmological 
constant problem, because gravity in 5D bulk is screened into 4D brane. But this model 
is disfavored by observations and suffers from a ghost problem. The ghost problem may 
be avoided by introducing a 6D theory or Galilean gravity. 

Inhomogeneous LTB model 

The inhomogeneous LTB model implies that current cosmic acceleration is not real at all. 
In other models, we have assumed FLRW metric, i.e. a homogeneous isotropic Universe. 
In the inhomogeneous LTB model we drop out the homogeneity, but there is a huge 
underdense void. If one measures a larger expansion rate (weaker gravity in underdence 
region) for lower redshift than higher redshift, it appears as cosmic acceleration. We may 
live in an underdense void which leads to this apparent cosmic acceleration. In large-scale 
structures we indeed observe inhomogeneous voids. But it is still challenging to satisfy 
CMB constraints. And this model implies that we are located at a special point (within 
15 Mpc from the center of void). 

The current status for probing dark energy is as follows. All current data provide 
strong evidence for an accelerated expansion of the Universe. SNe, CMB and BAO data 
provide independent evidence that the Universe contains 73% dark energy which only 
came to be dominant after the observed structure had formed. How can we explain this 
fact? It is an important issue in physics and cosmology today, what is the nature of dark 
energy and its time evolution if it exists. Is dark energy actually a fine-tuned vacuum 
energy and are we just lucky (as anthropic consideration)? Do quintessence, k-essence, 
coupled or unified models, modified gravity models or inhomogeneity models give better 
explanations? The k-essence model is one of the most self consistent models. The DGP 
model or the inhomogeneous LTB model may bring us a completely different view of the 
world in future. We do not have one promising model yet. Will the Universe recollapse 
or continue to expand with or without acceleration? The destiny of the Universe depends 
definitely on the nature of dark energy. 

There are many space- and ground-based observations in progress or being planned. 
They are aiming at selecting an appropriate dark energy model. It is still open, how 
we reduce systematic errors for observations. To get more reliable observational data we 
need huge volumes probed by the most recent deep redshift surveys as well. We have 
characterized dark energy and its effects on the expansion through the equation of state 
parameter w which is not what we actually measure. Thus it is convenient to introduce 
improved descriptions of dark energy which are better matched to observation. And we 
need theoretical models that make specific enough predictions to select out the better 
matching model or distinguish them from the ACDM model. For example, the chameleon 
mechanism provides very specific predictions for the effective Newtonian constant and the 
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Eotvos parameter in satellite measurement, as we discussed in section [5j 

From Einstein's biggest blunder to the Nobel Prize in Physics 2011: now we are ready 
to meet the most mysterious cosmological discoveries. 
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A Basics of Cosmology 
A.l Homogeneous FLRW model 

We will introduce here some cosmological quantities. The cosmological dynamics are given 
by solving the Einstein equation 

= 8nGT^, (183) 
where the Einstein tensor G^ is defined as 

G>2, = R pi/ — -g^ u R. (184) 

Assuming the Universe is spatially homogeneous and isotropic its metric takes the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) form 



ds 2 = g^ u dx^dx u = dt 2 — a 2 {t) 



<h ' 2 +r 2 dQ 2 



kr 2 



;i85) 



where dfl 2 = <i# 2 +sin 2 Odcf) 2 and the sign function of the curvature parameter k takes values 
+1, and -1 for positively curved, flat and negatively curved spatial section, respectively. 

We use the scale factor which is the relative size of the spatial sections as a function of 
time a(t) = R(t)/Ro, where the subscript represents the present time. The scale factor 
can be expressed in terms of the redshift z as a = normalized so that do = 1. 

On large scales matter in the Universe may be modeled as a perfect fluid which is char- 
acterized by energy density p, isotropic pressure p and four-velocity = -iA normalized 
so that u^u^ = 1. The energy-momentum tensor is given by 

Tp V = (p + p)upU v - pg^u- (186) 

The energy conservation equation {T a . a = 0, where ;Q denotes the covariant derivative 
with respect to a component) in an isotropic and homogeneous universe implies 

p = -3H(p + p), (187) 
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where the Hubble parameter is defined as 



H = -. (188) 
a 

Another convenient dimensionless parameter that characterizes the expansion is the de- 
celeration parameter: 

a 

» s -sp- (189) 

For the accelerating Universe q is defined to be negative. 

Substituting the metric and the energy momentum tensor into the Einstein equation, 
we get the first and second Friedmann equations: 



k 8ttG 
a 

and 



a 4ttG . . 

- = -— p + 3p. 191 
a 3 

For the relation between energy density p and pressure p, we define the equation of 
state as 

w = -. (192) 
P 

For example, non-relativistic cold matter (e.g. dust) is pressureless and corresponds to 
w — 0, radiation to w — |, and the cosmological constant to w = —1. Combining the 
equation of state and the second Friedmann equation, we get the condition —1 < w < — ^ 
for the observational status, namely the accelerated expansion of the Universe. In general 
relativity, the null dominant energy condition ensures that energy does not propagate 
outside the light cone, i.e., \p\ < \p\ which leads to — 1 < w < 1. However, we see that 
some models violate such conditions, phantom {w < —1) or superluminal sound speed 
(w > 1). 

For a flat universe (k = 0), solving the Eqs.(187) and (190) or (191), we express 
cosmological quantities in terms of w as 

P««- 3(1+w) , a(*)cxtTO, H= 2 (193) 

3(1 + w)t 

where w is assumed to be time independent and w ^ — 1. For w — — 1, 

p = const., a(t) oc e Ht , H = const. (194) 

The total energy density can be divided into four components (cold matter, radiation, 
dark energy and curvature density) 

Ptot = p m + Pr + Pa + Pk (195) 
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and each component decays as powers of the scale factor a: 



Ptot — Pmfl 



where p^ = — 3/c/(87rGa 2 ) and p ii0 = Pi(t ). Here we have assumed w — — 1 for dark 
energy component. One can see that pa becomes relatively important as time goes by. 
Now we define the present density parameters as 

n< - (197) 

Per (to) 

where the critical density p cr = j^S is the total energy density of matter or energy needed 



for the Universe to be spatially flat. It follows from Eqs. (190) and (196) that 

i = J2 n ^ ( 198 ) 

i 

where Qf. = —k/^aQH^) and f4 ~ from figure [T7T] 

In the Eq. ( 195[ ) we have introduced the ACDM model, which is impressively consistent 



with current observations. The ACDM model is a cold dark matter model with dark 
energy, which has been proposed to explain observational results such as anisotropics of 
the CMB, the large scale structure of galaxy clusters, accelerating expansion and chemical 
distribution of the Universe. 

The letter A stands for the dark energy component which occupies 72.8^Jg% of the 
total energy density of the Universe. The cold dark matter is non-relativistic and non- 
baryonic, and it does not interact with itself or other particles but interacts only through 
gravity. It is thought as making up 22.7±1.4% of our Universe. The remaining 4.56±0.16% 
is baryonic matter which builds visible planets, stars and galaxies. The current Hubble 
parameter is H = 70A~^\'\ km/s/Mpc (From WMAP+BAO+ifo data which are released 
at 2011 [152J). 

The ACDM model is built up on the homogeneous, isotropic and flat Universe. In 
the ACDM model, the Universe was hot and dense and had an early phase of expansion, 
when the light elements are produced via big bang nucleosynthesis and the radiation 
plasma dominates the total energy density of the Universe. Before the radiation dominated 
phase, the Universe had an earlier epoch of accelerated expansion, called inflation, when 
primordial density perturbations were seeded by quantum fluctuations, leaving an imprint 
on the CMB anisotropy and leading to the formation of structures due to gravitational 
instability on the present, isotropic and flat Universe. 

A. 2 Scalar field theory 

There are several dark energy models which are based on a scalar field such as quintessence 
and k-essence. In this subsection we introduce the basics of scalar field theory which we 
need for understanding such models. 
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A scalar field is a field which is invariant under any Lorentz transformation and has a 
scalar value at any point in space. For example, the temperature of a swimming pool is a 
scalar field, i.e., at each point we can measure a scalar value of temperature. The scalar 
field is known as spin-zero particles, and recently the Higgs particle which is a scalar field 
has been observed at LHC. Scalar fields are often introduced because of mathematical 
simplicity. To combine classical fields and general relativity we use the principle of the 
least action. The action S is the time integral of the Lagrangian L 

S[xi] = J L[ Xi (t),Xi(t)]dt (199) 

and can be rewritten by integrating the Lagrangian density £ over all space-time 

S[<pi] = [ C[^{x)]d A x, (200) 



where <fi(x) are classical fields including scalar and metric fields. In classical mechanics 
the Lagrangian is defined as the kinetic energy of the system minus its potential energy 



T-V. 



(201) 



For classical fields, the Lagrangian is very useful because it is invariant under some trans- 
formations that represent symmetries of the system. We can directly obtain the Euler- 
Lagrange equation by taking the least action : 

SS dL d(3L\ . 



(203) 



Sxi dxi dt \ dii 
The Einstein-Hilbert action describing gravity in general relativity is given by 

1 



S = 



16nG 



Ry^gd x, 



where g = det(g^ u ). Then the full action is the sum of the Einstein-Hilbert term and a 
matter (e.g. scalar fields) term, C M = — V(0), 



S = 



16nG 



R + C 



M 



gd x. 



(204) 



Note that \f—g is introduced in the integration measure to make C-^^g a scalar density, 
and gravity couples to matter minimally through this factor. Einstein gravity assumes 
this minimal coupling to ensure the equivalence principle. 

According to the action principle, the variation of an action with respect to g^ v is zero 



SgltfS 



1 5{^R) , 5(v^£ 



16nG 
1 



5R 



+ 



Mi 



R S(y/=g) 



(205) 



bg^d^x 

l S(y^gC M ) 



'-9 



Sg 



[IV 



Sg^V=9d 4 x. 
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As it is valid for any variation 5g^ v ', the square bracket has to vanish. Thus we obtain the 
following equation of motion 



SR 



+ 



R S 



9 



16nG- 



1 5(y/=gC 



Mi 



Sg pv y/^g Sg pv yj-g Sg pv 

We define the right hand side of the equation as an energy momentum tensor T M „ 

2 6{y/=gC M ) 



T, 



flV 



^ Sg 



ft i' 



(206) 



(207) 



From this equation of motion we derive the Einstein's field equation. To do that we need 
to know what 5R and SyJ—g are. We start with the Riemann tensor in order to calculate 
SR. 



dJT 



p 4- ^r p r A -+- r p ^r A 

Ha + 01 fi\ ua^ 1 fiX 0L va 



r p sr x 



(208) 



and the variation ST? is the difference of two connections, therefore it is itself a tensor 



[153J. Thus its covariant derivative is given by 



v„(ar' 



From this expression we can directly read off 

SR^ U = V,(ST P 

=► sr, u = v p (ary - v„((jr^ 

The Ricci scalar is defined by R = g^R^, then 

5R = R^SgT + g^SR^ 

= R pv Sg^ + VpfTSTP - V v g^SY 



(209) 
(210) 



(211) 



R^Sg^ 



The second term of the right hand side is a total derivative, and thus, according to Stoke's 
theorem, it can be expressed as boundary term by integral. As we assume the variation of 
the metric Sg^ vanishes at infinity, this term does not contribute to the variation of the 
action. Therefore, we get the following expression 



SR 



fig 



R 



■\iv 



(212) 



Now we consider Sy/^g. Since Sg = Sdet(g IIU ), according to Jacobi's formula^ it can 
be rewritten as gg pu Sg^ v . Then it follows 

1 . 1 
2 



-9 



-9 



--Sg 



-g(g" u Sg, 



flU J 



(213) 



3 The Jacobi's formula is a rule for differentiating a determinant. If a matrix A is invertible for all t, 
then f t detA(t) = detA(t)tr (A^)' 1 f t A{t)) . 
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Since S^g^g^) = 0, it follows that g^Sg^ = —g^Sg^ and 



/IV 



(214) 



The left hand side of the equation of motion (206) finally becomes R^ u — \g^ v R so that 
the equation of motion takes the form of Einstein's field equation 



Rfj,v — -^9»vR = SttGT^. 
We come back to the energy momentum tensor 

2 5{y/=gC M ) 



T, 



-g 8g 



The variation with respect to 5g^ u splits into two parts, ^J= S ^J 1 Cm and ^=|p6- 



(215) 



(216) 



-g- 



From the variation of the determinant (|214|), we get 

-g^u^M- 



2 5^ 



-g V 



c 



A I 



(217) 



Now we assume a scalar field as the matter sector, i.e., Cm — ^d fl (f)d fl (j) ~ ^(0)> then 



2^=2- 6 



g ap d a 



- 2- 



5 



-v{4>). 



Hence 



& a ti&^v and V(4>) being independent of the metric g^ v leads to 
n 8C M 



(218) 



Sg 



fJbV 



6ffid a <l>dp<l> 



(219) 



0, 



Now we have the energy momentum tensor in term of the scalar field 

"1 



9i*u 



-d a <pd a <p - V{<f>) 



or equally 



T\ = &><}>d v <t> - 



-d a <f>d a <f>-V(<f>) 



5\ 



(220) 



(221) 
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For the (0,0) component we get T° — p — \4> 2 +V(<f>). For the (i,i) component we 
get —\T % i = p = \4> 2 — V{4>). If we take the FLRW metric and assume the Universe is 
homogeneous, we obtain following Friedmann equations 



B'+^ = ^(^ + V(*)), (222) 



a 8nG 



¥ - V{$)) ■ (223) 



a 3 

Let us check again the Lagrange density C. The full Lagrange density including the 
gravity term and the scalar field term is 

C ^~ 9 = ~iLg R ^~ 9 + l^tV^g - V{ct>)^r g . (224) 
The functional derivative of C with respect to <f> is zero according to Euler lagrange equation 



5S[g, 



v d{d v <t>) 

-V'i&yf^j - {d v d v 4>)^-g - V<l>{d v y/=g) = 0, 



(225) 



where the prime denotes 

We take again the FLRW metric, and assume the homogeneous Universe to drop 
spatial gradients. For a diagonal matrix, the determinant is the product of the diagonal 
components 

a 6 r 4 sin 2 8 

9 = 



1 — kr 2 

Because of spatial homogeneity d l (p vanishes, and we obtain 





■V'WV^g - ^V^g - <t>{d t y/=g). 



Since dt^f^g = = —y/—g, we obtain 



2 a 

O = -U'(0)-0-0-, 

(X 

which yields the scalar field equation 

4> + 3H<j) + V'((f)) = 0. 



(226) 



(227) 



(228) 



(229) 



In the same way, we can derive the equations of motion from an action for each model. 
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